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Abstract. Weakly interacting dark matter particles can be pair-produced at colliders and detected through 
signatures featuring missing energy in association with either QCD/EW radiation or heavy quarks. In 
order to constrain the mass and the couplings to standard model particles, accurate and precise pre¬ 
dictions for production cross sections and distributions are of prime importance. In this work, we con¬ 
sider various simplified models with s-channel mediators. We implement such models in the Feyn- 
Rules/MadGraph5_aMC@NLO framework, which allows to include higher-order QCD corrections in 
realistic simulations and to study their effect systematically. As a first phenomenological application, we 
present predictions for dark matter production in association with jets and with a top-quark pair at the 
LHC, at next-to-leading order accuracy in QCD, including matching/merging to parton showers. Our study 
shows that higher-order QCD corrections to dark matter production via s-channel mediators have a sig¬ 
nificant impact not only on total production rates, but also on shapes of distributions. We also show that 
the inclusion of next-to-leading order effects results in a sizeable reduction of the theoretical uncertainties. 


1 Introduction 

Various astrophysical and cosmological observations pro¬ 
vide strong hints for the existence of dark matter (DM). 
Yet, very little is known about the nature of DM, and con¬ 
straints on models from various direct/indirect detection 
experiments and cosmology still allow for a wide range of 
DM masses and couplings to the Standard Model (SM) 
particles. Various types of DAI searches are sensitive to 
different regions of the DM model parameter space |l|2|. 
In order to maximise the chances for discovering - or at 
least excluding - DM models, it is hence imperative to per¬ 
form both astrophysical and collider searches for DM. The 
most promising way to detect signals of weakly interact¬ 
ing DM particles at the LHC is through their associated 
production with jets, EW bosons and heavy quarks, lead¬ 
ing to signatures with missing transverse energy (MET). 
Searches for DM have been performed at LHC Run I (see 
e.g. Refs. 3|[4]) and are one of the central goals of LHC 
Run II 5 |. 

While the complementarity of different DM searches 
is a powerful tool, it is intrinsically limited in that the 
comparison of collider and other results introduces some 
degree of model dependence. As the nature of DM is still 
unknown, there exist a myriad of DM models and mecha¬ 
nisms to be explored, spanning a wide range of complexity 
and ambition. In this context, it is of utmost importance 


to follow an approach where the model dependence is lim¬ 
ited, while the salient features that can provide a useful 
characterisation of possible signals are kept intact, i.e., 
the “simplified models” approach [6 : . 


Many simplified models for DM have been proposed 
in the past (see Refs. [5,7,8] and references therein). In 
their simplest form, these models assume DM to be a sin¬ 
gle massive particle which interacts weakly with the SM 
particles. The interaction of DM with the SM can be me¬ 
diated by a new field, which we dub the mediator. When 
the mediator is heavier than the energy scales the exper¬ 
iment can probe, the interaction becomes point-like and 
the Lagrangian can be organised in terms of a tower of 
higher-dimensional operators in the framework of an effec¬ 
tive field theory (EFT). However, when the experimentally 
accessible energy becomes comparable or higher than the 
mediator mass, on-shell effects become important and a 
properly defined quantum field theory featuring the medi¬ 
ator state in the spectrum is needed [9 14. The LHC can 
explore a large range of DM and mediator masses, as well 
as coupling strengths and possible combinations of DM 
and mediator spins. Collider results within the framework 
of simplified models can then be combined with direct 
searches for the mediators, e.g. in Drell Yan processes, as 
well as with cosmological and astrophysical constraints on 
DM. 
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LHC searches for DM will rely on precision calculations 
to impose the most stringent bounds on DM models, and 
hopefully characterise possible signals. Higher-order cor¬ 
rections in QCD to DM production processes at the LHC 
are hence vital for extracting precise information about 
DM from the LHC results. 

Previous analyses 


15-19 studied next-to-leading or¬ 


der (NLO) QCD corrections for DM production in MET 
+ j/ r y/W in the case of EFT, i.e., in the limit of heavy 
mediators^ In this article, we consider simplified mod¬ 
els with s-channel mediators. We analyse the impact of 
the higher-order corrections to mono-jet signals in various 
benchmark scenarios with spin-1 (vector and axial-vector) 
mediators, calculate DM production cross sections (both 
total and differential) at NLO accuracy with up to two 
jets, and also merge the corresponding samples via the 
FxFx procedure [20]. To our knowledge, such accuracy 
has so far not been achieved in the context of DM simu¬ 
lations/predictions for the LHC Run II. 

In addition, we consider tt associated production and 
compute NLO cross sections and distributions for spin- 
0 (scalar and pseudo-scalar) mediators in representative 
cases, including those with a light mediator. Predictions 
for this class of processes at NLO in QCD also represent 
a novelty in the context of DM computations for LHC 
Run II. 

The first goal of this work is to illustrate the feasibility 
of having a fully general i mple mentation of DM simplified 
models in the FeynRules [21|/MadGraph5_aMC@NLO 
[22] (MG5aMC henceforth) framework, accurate up to 
NLO in QCD. To this aim we start with the simplest 
(yet non-trivial) case of s-channel mediators (colour sin¬ 
glet, spin 0 and 1 bosons) coupling to DM and quarks. 
We assume that DM is a Dirac fermion for concreteness, 
yet our implementation is not limited by the choice of 
DM spin or nature (Dirac or Majorana). We show how 
predictions and event generation for this class of models 
can be achieved at NLO QCD accuracy, in a fully auto¬ 
matic way, for a wide set of observable/final states, while 
also employing the most recent matching/merging multi¬ 


jet techniques 20 


The second goal of this work is of phenomenological 
nature, i.e. to investigate the impact of the NLO correc¬ 
tions on the current and future searches for DM at the 
LHC. We consider two examples, among several promis¬ 
ing final state signatures: 


pp —> XX + jets 


( 1 ) 


for a spin-1 mediator model, and 
pp -A XX + tt 


( 2 ) 


for a spin-0 mediator model, where X is a DM particle. We 
do not only calculate NLO QCD corrections to the overall 
production rates, but also study in detail the effects of 

1 A first discussion of the impact of NLO QCD corrections 
to DM production in a specific simplified model can be found 


in 17 and more recently in [5 18 . 


higher-order corrections on the differential distributions 
of relevant observables. 

Our simulation set-up is based on the FeynRules/ 
MG5aMC framework. The FeynRules package provides 
the relevant Feynman rules starting from any local La- 
21123 ) 


, as well as the UV/ R 2 counterterms [25] 
necessary for the NLO computations via Nloct 26]/Fey- 
nArts 27 . Our simplified DM model files are publicly 

With these 


grangian 


24 


available at the FeynRules repository 28 


ingredients, which are only based on the model and are not 
process dependent, MG5 aMC c omputes tree-level am¬ 
plitudes, loop-amplitudes 


for a desired process 
phase space 


29 

31 

we 

as 


and subtraction terms 
as their integration over 


32 . Event generation is obtained by match¬ 


ing short-distance events to the parton shower employ¬ 
ing the MC@NLO method 33 , which is implemented for 
Pythia6 [34], Pythia8 1351, Herwig 6 [36] and Her- 


WIG-| —b 371. We note that our DM model files can be 
exported not only to event generators but also to tools for 
DM relic abundance as well as DM direct and indirect de- 

38|3 9 


tections such as MicrOMEGAs 


and MadDM [40, 
41 , allowing for more comprehensive DM studies. 

The paper is organised as follows. In Sect. [2] we intro¬ 
duce simplified models for DM and specify the relevant in¬ 
teractions and model parameters of the implementation. 
We discuss the impact of the NLO QCD corrections on 
DM pair production with jets through spin-1 mediators 
in Sect. [3] which includes a discussion of inclusive cross 
sections, differential distributions as well as a discussion 
of the impact of parton showers and the NLO merging of 
events with different jet multiplicities. NLO QCD predic¬ 
tions for DM production in association with a top-quark 
pair through spin-0 mediators are discussed in Sect. [4] We 
provide our conclusions and an outlook in Sect. [5] 


2 Simplified dark matter models: 
the s-channel mediator case 

We start by defining the particle content and the inter¬ 
actions of the simplified model, which we dub DMsimp. 
We assume that DM is described by a single, massive and 
weakly interacting particle, that communicates with the 
SM through the exchange of a mediator. For simplicity, 
we assume that the mediator is not part of the SM^The 
first very general classification stems from the class of ver¬ 
tices that characterise the model: Lagrangians featuring 
DM-DM-mediator and SM-SM-mediator type interac¬ 
tions identify models with an s-channel mediator, while 
Lagrangians characterised by DM-SM-mediator interac¬ 
tions define a t-channel mediator. The former interactions, 
for example, arise in models featuring an extra scalar or 
Z' which couples to a pair of DM particles, while the lat¬ 
ter is common in supersymmetric models. From the point 
of view of QCD corrections the two classes are very dif¬ 
ferent, as an s-channel mediator is necessarily a colour 
singlet, while a t-channel mediator can be either neutral 

2 So-called portal models where DM interacts with the Higgs 
or the Z /7 bosons can be included in our framework as well. 
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or coloured. In this work we focus on the s-channel mod¬ 
els, leaving the implementation and validation of f-channel 
models to forthcoming work. The s- and f-channel classes 
can be further organised by the quantum numbers of the 
DM particle and the mediator. To start with, we focus on 
the case of Dirac DM with spin-1 or spin-0 mediators cou¬ 
pling to the matter fields of the SM. Changing the spin 
or the nature of the fermion (Dirac or Majorana) of the 
DM particle or including a coupling of the mediator to the 
SM bosons is straightforward [42]. On the other hand, ex¬ 
tending our analysis to spin-2 mediators, while feasible in 
principle, entails dedicated validation work, as such mod¬ 
els are, in general, not renormalisable. We defer such an 
extension to the future. 


2.1 Spin-1 mediator model 

In the framework of our simplified model, the interaction 
Lagrangian of a spin-1 mediator (Yi) with a Dirac fermion 
DM (Xd) is given by 

C^ D = x d1 ^ 9 v Xd + gi Dl5 )X D Y? , (3) 

and with quarks by 

£ SM = J2 

ij 

+ UiJvidua +9uij TsK 

where d and u denote down- and up-type quarks, respec¬ 
tively, (i,j= 1,2,3) are flavour indices, and g v ' A are the 
vector/axial-vector couplings of DM and quarks. Note that 
we adopt this notation according to the actual imple¬ 
mentation in FeynRules. The model file, including an 
alternative choice for the spin of DM particle (complex 
scalar Xq), can be downloaded at the FeynRules repos¬ 
itory [28] . 

The pure vector and pure axial-vector mediator sce¬ 
narios are given by setting the parameters in the La- 


grangians 

§ 

and Q to 


9x d = 

9x 

and g% D = 0 

(5) 

9uu = 

9d u 

= g SM and g A . = gf.. = 0 

(6) 

and 




9x d = 

0 

and gj Co = g x 

(7) 

9u« = 

9du 

= 0 and g A . = g A _. = g SM , 

(8) 


respectively, where we assume quark couplings to the me¬ 
diator to be flavour universal and set all flavour off-diagonal 
couplings to zero. With this simplification of a single uni¬ 
versal coupling for the SM-Yi interactions, the model has 
only four independent parameters, i.e. two couplings and 
two masses: 


1 1 • 


(4) 


(9) 


We note that the mediator width is calculated from the 
above parameters. 

Finding a signal of DM in this parameter space (or 
to constrain these parameters) is the primary goal of the 
DM searches at the LHC Run II [5], and the most impor¬ 
tant signature in this model is mono-jet plus MET. The 
di-jet final state via the Y) Drell Yan process can be an 
important complementary channel. 


2.2 Spin-0 mediator model 


Similarly, in the case of a spin-0 mediator (Y 0 ) interact¬ 
ing with the Dirac fermion DM and the SM particles, we 
define the interaction part of the Lagrangians as 

Cl? D = X d { 9 s Xd + igx D ls)X D Fo , (10) 

and 


rY 0 

■‘■'SM 


H [dr^ial, +i9i j l5)d J 


l ^2 [9u 


igLis ) u j 


Yo, 


( 11 ) 


where g s ^ p are the scalar/pseudo-scalar couplings of DM 
and quarks. Assuming a UV complete description of the 
scalar theory with the couplings of the mediator to the 
SM particles proportional to the particle masses, we nor¬ 
malise these couplings to the SM Yukawa couplings, y{ t = 


y/2m.f/v, and set all flavour off-diagonal couplings to zero. 
This implies that, in a five-flavour scheme with massless 
bottom quarks, only top quarks are relevant for DM pro¬ 
duction in this model. Extension to a four-flavour scheme 
with massive bottom quarks is possible. The model file for 
the spin-0 mediator case, including other choices for the 
spin of the DM particle (real scalar Xr and complex scalar 
Xc), is also available at the FeynRules repository 


28 


The pure scalar and pure pseudo-scalar mediator sce¬ 
narios are given by setting the parameters in the La¬ 
grangians (10) and {H]> to 


9x d 

= gx 

and 

9x n 

= 0 

(12) 

9 S 

Yu 33 

= 9SM and g p 33 

= 0 

(13) 

and 






9Xn 

= 0 

and 

9x D = 

9x 

(14) 

9 s 

y« 33 

= 0 

and 

P . 

9u 33 — 

9sm , 

(15) 


respectively. All the other 9 u[f and g S J P parameters are 
irrelevant. Similar to the spin-1 case, the model has only 
four independent parameters as in ([9]). 

In the spin-0 mediator model with Yukawa-type cou¬ 
plings, the most relevant tree-level process at the LHC is 
DM pair production associated with a top-quark pair. On 
the other hand similarly to Higgs production, at one loop, 
gluon fusion can give rise to MET + jets signatures which 


{ 9 sm , gx, m x , m Y } ■ 
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are in general phenomenologically relevant. For the heavy 
mediator case, the four-top final state can be also relevant. 

At this stage, we do not see compelling reasons to in¬ 
troduce couplings of the mediator to leptons, even though 
it is straightforward to do. We do not include effective 
gluon-gluon-Fo interactions either, for several reasons. The 
first is that this operator is higher dimensional (dim=5) 
and therefore might lead to unitarity-violating effects that 
need to be studied on a model and benchmark basis. The 
second is that a simplified model assumes no other new 
physics particle beyond the dark matter particle X and 
mediator Y at the weak scale. If such new particle ef¬ 
fects decouple with their mass, the main contribution to 
the gluon-gluon-lb coupling would then be due to loop 
of SM quarks. Depending on the masses involved (that 
of the DM, the mediator and the quarks contributing in 
the loop) and the momentum transfer of the process, such 
interactions might be considered point-like. This is, how¬ 
ever, very much model and process dependent. To be accu¬ 
rate one should first calculate the loop-induced processes 
exactly and study the range of applicability of such an ef¬ 
fective interaction by explicit comparison. This is possible 
in MG5 aMC 43 and has been considered with the same 


DM implementation as presented here in Ref. 44 . Other 
studies of the loop-induced process for mono-jet + MET 
can be found in Refs. 45 ■ 48]. We also note that couplings 
of the mediator to the SM gauge bosons can be introduced 


easily 42 


3 Dark matter production with jets 

In this section, we consider a spin-1 mediator scenario and 
discuss the impact of the NLO-QCD corrections on DM 
pair production with jets, be., 

pp^XX + j(j). (16) 

In MG5 aMC the code and events for the above process 
can be automatically generated by issuing the following 
commands: 

> import model DMsimp_s_spinl 

> generate p p > xd xd~ j [QCD] 

> add process p p > xd xd~ j j [QCD] 

> output 

> launch 

We have checked that our model can reproduce the SM 
predictions for pp —> Zj(j) —> t + t~ j(j) by adjusting the 
corresponding coupling and mass parameters. 

To illustrate the effect of the higher-order corrections, 
we consider pure vector, Eqs. © and ([6]), or pure axial- 
vector, Eqs. ([7]) and Q, couplings with a simplified flavour 
structure. We take 

9x = 1 and g SM = 0.25 (17) 

as our benchmark for the spin-1 mediator scenario. 

We assume that the mediator can only decay into the 
DM particle and the SM quarks (if kinematically allowed) 


10 3 

10 2 


-Q 

a, io° 

o , 
j 10 1 

2; 

b ID ' 2 


io~ 3 


10“ 



f = 50 GeV 


pp^XX+j 
at the LHC13 


: vector mediator 
r(i?X’ 0 sm) = (1’ 0-25) 


m x = 500 GeV 




m Y [GeV 


Fig. 1. Summary plot of NLO cross sections and corresponding 
K factors in Table |T] 


through the interactions specified in Eqs. © and 0, 
so that the mediator width, TV, is determined by the 
particle masses and the couplings <jx and gsM- In our 
framework, the width is automatically computed by using 
the Mad Width module 49 for each parameter point. 


The above benchmark coupling strength in © leads to 
Fy/my ~ 0.05 for my > 2mx and TV /my ~ 0.025 for 
uiy < 2 mx, both for the vector and axial-vector cases. 
Note that, if we take g$M = 1 (0-5) with g x = 1, the 
Y\ width becomes very large as TV /my ~ 0.5 (0.15) for 
my > 2 mx- 

We provide LO and NLO QCD predictions for pp —► 
XX +j(j) at the center-of-mass energy y/s = 13 TeV. The 
central value po for the renornralisation (pr) and factori¬ 
sation ( pp ) scales is set to Hr/ 2, where Hr is the sum 
of the transverse momenta of all jets in the event and the 
missing transverse energy. The scale uncertainty is esti¬ 
mated by varying the scales pr and pp, independently, by 
a factor two around po. We adopt the five-flavour scheme 
and the LO and NLO NNPDF2.3 set [50] through the 
LHAPDF interface 51 , with the corresponding values of 
a^°(Mz) = 0.130 and a^ LO (Mz) = 0.118, for the LO and 
NLO predictions, respectively. The PDF uncertainties are 
computed automatically [52], following the prescription 


summarised in 53 


Where relevant, we apply a parton shower to the events 
using Pythia8 35 . We then define jets using the anti -kp 

with the 


55 


algorithm [54j as implemented in Fast Jet 
jet cone radius R = 0.4, where we require Pt(JT> 30 GeV 
and | r](j) | < 4.5 for all jets in the event. 


3.1 Total cross sections 

In Table [T] we present LO and NLO cross sections (in 
pb) for DM pair production in association with a jet at 
fixed order (FO) in perturbation theory. We show results 
for a pure vector mediator by fixing the parameters as 
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K,m x ) [GeV] 



MET > 150 GeV 

vector 

MET > 300 GeV 

MET > 500 GeV 

10 

undecayed 

olo [pb] 

ONLO [pb] 

K factor 

2.923 x 10 2 ±?° 9 7 ±i. 6% 
5.093 x 10 2 ^8°2 3 ±o.5% 
1.74 

1.734 x 10 1 ±11.9 ±1.1 % 
2.689 x 10 1 ±9 °i 4 ±O.6% 
1.55 

1.695 x 10°±iI; 4 o ±i.8% 
2.433 x 10° ±io!o ±1-1 % 
1.44 

(10,1) 

my > 2 mx 

olo [pb] 

CTNLO [pb] 

K factor 

1.605 x 10 2 ±?° 9 7 ±i. 6% 
2.818 x 10 2 ±^±0.5% 
1.76 

0.978 x 10 1 ±12.0 ±i.i % 
1.517 x 10 1 ±8° 9 °±o.6% 
1.55 

0.970 x 10° ±11;? ±2.0% 
1.345 x 10°±?° 6 5 ±i.i% 
1.39 

(10,50) 

my < 2mx 

olo [pb] 

ONLO [pb] 

K factor 

2.434 x 10°tlo.?±i.5% 
3.198 X 10° ±54 ±0.5 % 

1.31 

2.843 x 10" 1 ±{2.5 ±12% 
3.485 x 10“ 1 ±6.3 ±o.7% 
1.23 

3.786 x 10 —2 ±{±±±2.4% 
4.325 x 10“ 2 ±±|±i.3% 
1.14 

(10,500) 

my < 2 mx 

olo [pb] 

ONLO [pb] 

K factor 

6.968 x 10~ 3 ±1±?±4.3% 
7.698 X 10" 3 ±6.4 ±2.2% 
1.10 

2.314 x 10 —3 ±15.0 ±4.6 % 
2.385 x 10“ 3 ±6.9±2.3% 
1.03 

7.317 x 10 —4 ±16.1 ±5.6 % 
6.800 x 10“ 4 ±?;?±2.6% 
0.93 

100 

undecayed 

olo [pb] 

ONLO [pb] 

K factor 

2.148 x 10 2 ±9°3 6 ±i. 5% 
3.011 x 10 2 ±5.9 ±0.5% 

1.40 

1.616 x 10 1 ±!a:<Ul.o% 
2.121 x 10 1 ±7.1 ±0.6 % 

1.31 

1.644 x 10° ±?I;? ±1.9% 
1.955 x 10° ±|;± ±1.2% 

1.19 

(100,1) 

my > 2 mx 

olo [pb] 

ONLO [pb] 

K factor 

1.100 x 10 2 ±9°3 6 ±i. 5% 
1.530 x 10 2 ±5* ±0.5% 

1.39 

0.822 x 10 1 ±12.0 ±i.i % 
1.100 x 10 1 ±7.2 ±0.6% 

1.34 

0.862 x 10° ±11? ±1.9% 
1.059 x 10° ±1? ±1.2% 

1.23 

(95, 50) 

my < 2m x 

olo [pb] 

ONLO [pb] 

K factor 

1.117 X 10 1 ±9 1 6°±1.5% 
1.512 X 10 1 ±5.5 ±0.5% 

1.35 

0.988 x 10° ±±±2 ±1.1% 
1.281 x 10° ±2;| ±0.6% 

1.30 

1.140 x 10 —1 ±H"_2 ±2.0 % 
1.325 x 10 —1 ±7.6 ±1.2 % 
1.16 

(100,500) 

my < 2 mx 

olo [pb] 

CTNLO [pb] 

K factor 

7.043 x 10" 3 ±1I;o±4.3% 
7.804 X 10 _3 ±6.4 ±2.2 % 
1.11 

2.329 x 10 —3 ±15.0 ±4.6 % 
2.411 X 10“ 3 ±6.8±2.3% 
1.04 

7.395 x 10 -4 ±?6.i ±5.6% 
6.908 x 10~ 4 ±?;?±2.6% 
0.93 

1000 

undecayed 

u"lo [pb] 

UNLO [pb] 

K factor 

2.248 x 10° ±13.2 ±3.2% 
2.601 x 10° ±6.o ±1-7% 

1.16 

6.865 x 10 —1 ±14.3 ±3.3 % 
7.393 x 10" 1 ±6.4 ±1.8% 
1.08 

1.979 x 10-^1^ ±4-1% 
1.909 x 10^ 1 ±6;g±2.i% 
0.96 

(1000,1) 

my > 2 mx 

u"lo [pb] 

ONLO [pb] 

K factor 

1.093 x 10° ±13.3 ±3.1 % 
1.215 X 10° ±5.5 ±1.7% 

1.11 

3.278 x 10" 1 ±11.4 ±3.3% 
3.399 x 10“ 1 ±6.o±i.7% 
1.04 

9.182 x 10“ 2 ±?5;6±4.i% 
8.743 x 10“ 2 ±?;s±2.o% 
0.95 

(1000,50) 

my > 2 mx 

<r lo [pb] 

UNLO [pb] 

K factor 

1.094 x 10° ±13.3 ±3.1 % 
1.221 X 10° ±5.6 ±1.7% 

1.12 

3.268 x 10 _1 ±}4;4±3.3% 
3.416 x 10" 1 ±6.0±1.7% 
1.05 

9.137 x 10“ 2 ±?5;6±4.i% 
8.807 x 10“ 2 ±?;6±2.o% 
0.96 

(995, 500) 

my < 2m x 

olo [pb] 

UNLO [pb] 

K factor 

2.169 x 10" 1 ±£3 ±3.4% 
2.497 x 10" 1 ±6.2 ±1-8% 

1.15 

6.777 x 10“ 2 ±}4.4±3.6% 
7.223 x 10“ 2 ±6.6 ±i.9% 
1.07 

1.981 x 10" 2 ±15.6 ±4.4% 
1.914 x 10“ 2 ±6;g±2.i% 
0.97 

(10000,1) 

my ^$> v § 

olo [pb] 

UNLO [pb] 

K factor 

8.487 x 10" 6 ±l4;3±4.3% 
8.835 x 10 —6 ±5.1 ±2.5 % 
1.04 

2.666 x 10 —6 ±15.7 ±5.5 % 
2.579 x 10 _3 ±5.5 ±3 % 

0.97 

8.238 x 10 —7 ±17.0 ±7.3 % 
7.148 x 10“ 7 ±lo±4.4% 
0.87 


Table 1 . LO and NLO cross sections and corresponding K factors for DM pair production in association with a jet for the 
vector mediator scenario at the 13-TeV LHC, where different MET cuts are imposed. The uncertainties represent the scale and 
PDF uncertainties in per cent, respectively. We show several benchmark model points for the mediator and DM masses with 
the coupling parameters gx = 1 and gsM = 0.25. 
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Fig. 2. MET distributions at FO (N)LO accuracy for pp —> XX + j at the 13-TeV LHC for (my,mx) = (100,1) and 
(1000,50) GeV, where we assume a pure vector mediator and Dirac DM. The middle and bottom panels show the differential 
scale uncertainties and I\ factors, respectively. 


in Eqs. ([5]), ^ and ( fl7| ). We cover various benchmark 
points suggested by theWTL AS/CMS DM forum |[5] in the 
my — tux plane, representing four different cases: on-slrell 
(my > 2m x) and off-shell (my < 2m x) production of the 
mediator, in the threshold regime {my < 2m x) and in the 
EFT limit (my \/S). We also present scale and PDF 
uncertainties in % as well as I\ factors which we define as 
the ratio of the central values of the NLO and LO cross 
sections. We compute the table entries with different MET 
cuts: 150, 300, and 500 GeV. For convenience, we also 
show a graphical summary of our results in Fig. [I] As a 
reference, the cross sections for pp —» Y\ +j are also shown, 
where the vector mediator is produced on-shell and 
does not decay. For my > 2mx, the mono-jet rate is given 
by cr{pp —> XX + j) ~ cr{pp —> Y\ + j) x B(Y\ —> XX) in 
the narrow width approximation. 

The production rate strongly depends on the both 
masses as well as on the kinematic cuts, and varies by 
orders of magnitude in the parameter scan. On the other 
hand, the K factors, i.e. higher-order effects, are not so 
sensitive to the mass spectra; K ~ 1.1 for the heavy- 
mediator and/or heavy-DM cases, while K ~ 1.3 — 1.4 
for the ~ 100 GeV mediator with light DM, assuming the 
MET > 150 GeV cut. We find that in the case of a rela¬ 
tively light mediator with a very light DM, ( my,mx ) = 
(10,1) GeV, the K factor can reach a value as large as 
1 . 8 . 

Different benchmark points probe different Bjorken- 
x regions of the parton distribution functions. As heavy 


mediators/DM are produced from very high-a; partons, the 
dominant contribution comes from the qq initial state, as 
the gluon PDF is sub-dominant in the high-a; region. For 
light mediators with light DM, on the other hand, a large 
contribution arises from the qg initial state. For instance, 
we find that the ratio of production cross sections via uu 
and ug initial states, a{uu)/a{ug), is 1.4 in the case of 
{my,mx) = (1000,50) GeV while 0.2 for ( my,mx ) = 
(100,1) GeV, at LO. 

As expected, most of the results at NLO accuracy 
display significantly smaller scale uncertainties compared 
to the LO calculations. An exception is provided by the 
( my,mx) = (10,1) GeV case, which we discuss in detail 
at the end of this subsection. The PDF uncertainties are 
sub-leading in both the LO and NLO results and reduced 
by going from LO to NLO. Furthermore, the scale and 
PDF uncertainties increase when the mass scale of the 
mediator and/or DM increases. 

The higher MET cut leads to smaller I\ factors and to 
larger scale and PDF uncertainties, which one can clearly 
see in the MET distributions in Fig. [2] 

In Table [2j we present the pure axial-vector med iato r 
case by fixing the parameters as in Eqs. Q, (|8]) and (17). 
The resulting cross sections are very similar compared to 
the pure vector case for my > 2mx, while in the off-shell 
regime, we find that the cross sections are suppressed com¬ 
pared to the production via pure vector mediators. In the 
off-shell situation the DM pair is produced at threshold. 
A pair of (Dirac) DM originating from a decay of a spin-1 
mediator will be in a 2 S+ 1 Lj state with J = 1. If the 
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(m Y ,m x ) [GeV] 



MET > 150 GeV 

axial-vector 

MET > 300 GeV 

MET > 500 GeV 

100 

undecayed 

a LO [pb] 
ONLO [pb] 

K factor 

2.130 x 10 2 ±g° 3 6 ±i.6% 
3.063 x 10 2 ±6.i ±0.5% 
1.44 

1.573 x 10 1 ± 12.0 ± 1.1 % 
2.153 x 10 1 ±7)1 ±0.6% 

1.37 

1.633 x 10° ±±I.o ±1.9% 
2.055 x 10° ± 8.3 ± 1 . 6 % 

1.26 

(100,1) 

my > 2 mx 

CLO [pb] 
ONLO [pb] 

K factor 

1.101 x 10 2 ± 9 ° 3 6 ±i. 6 % 
1.549 x 10 2 I®.® ±0.5% 

1.41 

0.825 x 10 1 ± 12.1 ± 1.1 % 
1.127 x 10 1 ±7.2 ± 0.6 % 

1.37 

0.854 x 10°±il;t±2% 

1.063 x 10° ±l;a±i.2% 

1.24 

(95,50) 

my <2mx 

cmo [pb] 

ONLO [pb] 
K factor 

3.070 x 10° ±iao ±1-5% 
4.093 x 10° ± 5 ° ± 0 . 5 % 
1.33 

3.359 x 10 _1 ± 12.4 ± 1.2 % 
4.302 x 10 - 1 ±e: 9 ±o. 7 % 
1.28 

4.457 x 10 -2 ±il '.3 ± 1 . 8 % 
5.079 x 10 ~ 2 ±?; 4 ±i. 3 % 
1.14 

(100,500) 

my < 2 mx 

a ‘LO [pb] 
ONLO [pb] 
K factor 

2.298 x 10 - 3 ±]l 5 ± 5 % 
2.502 x 10 " 3 ± 6.8 ± 2 - 5 % 
1.09 

7.839 x 10- 4 ±i 9 5 ±±5.3% 
7.972 x 10 “ 4 ± 7 ; 3 ± 2 . 6 % 

1.02 

2.558 x 10 - 4 ±?e: 5 ± 6 . 3 % 
2.383 x 10 " 4 ±?; 5 ± 3 .o% 
0.93 


Table 2. Same as Table 1 but for the axial-vector mediator scenario. 



0 50 100 150 200 


MET [GeV] 



Fig. 3. MET distribution at (N)LO+PS accuracy for pp —> XX + j at the 13-TeV LHC for (my,m.v) = (10,1) GeV (left) 
and ( m,Y,mx ) = (100,1) GeV (right). The lower panels provide information on the differential scale uncertainty and K factor. 


coupling is vector-like the DM pair can be in a 3 Si state, 
while if it is axial-like it will be in a 3 Pi state, i.e. sup¬ 
pressed at threshold. The similar argument holds in case 
of gg —¥ Y 0 + tt case, as we show in Sect, [fj 

The NLO effects are very similar to the vector me¬ 
diator scenario for all the mass combinations as well as 
the MET cuts. Although we do not show the mixed sce¬ 
nario of vector and axial-vector, one can easily compute 
such scenarios by changing the coupling parameters in our 
simplified model. 


The parameter point (my,m.y) = (10,1) GeV war¬ 
rants special attention, as it illustrates a case of large NLO 
corrections (so-called “giant K factors” [56]), which might 
arise in the limit where p? T m Y ,mx- The giant K fac¬ 
tors in the pp Y\ + j process occur due to the opening 
of the pp —► Y-\ + jj channel at NLO. This process can 
lead to a di-jet event topology with a soft, collinear emis¬ 
sion of Y\. In the regime of p J T m Y , the Y\ emission 
behaves similar to an emission of a massless gauge bo¬ 
son, where the diagrams with dijet topologies contribute 
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factors of a 2 s g 2 x \og? (pUp / my), and hence NLO I\ factors 
which scale as ~ a s \og 2 (p J T /my). Similar features com¬ 
monly appear in calculations of SM W/Z+ jets processes 
at high jet py 


56 


Topologies leading to giant K factors are naturally 
suppressed in the case of DM production by the cut on 
MET. This restricts the calculation to regions of phase 
space which are insensitive to the soft and collinear double¬ 
logs of di-jet event topologies with a soft Y\ emission. Fig¬ 
ure [3] (left) illustrates the effect in case of my = 10 GeV 
and nix = 1 GeV. The region of low missing energy dis¬ 
plays a two to three orders of magnitude difference in 
rate between the LO and NLO calculations, whereas we 
see that above MET > 50 GeV, the K factor is drasti¬ 
cally reduced. On the other hand, we see that already for 
(my, nix) = (100,1) GeV in Fig. [ 3 ] (right), such logarith¬ 
mic enhancements are only very weak. 


3.2 Differential distributions 

We proceed with the discussion of the features of the dif¬ 
ferential distributions relevant for DM studies. We begin 
with Fig. [f] which shows the MET distributions at LO and 
NLO for four benchmark points of the simplified model, 
assuming a pure vector mediator and Dirac fermion DM. 
As seen in the total rates, the NLO effects in the distri¬ 
butions do not depend on the mass relation between the 
mediator and the DM, i.e. on-shell or off-shell, but do 
depend on the energy scale of the final state. In the top 
panels, the energy scale is 0(100) GeV for my or 2 mx- 
We find that the two benchmark points display striking 
similarities in the shape of the MET distributions, while 
the rate of the latter is suppressed due to off-shell Yi pro¬ 
duction. The largest effects of NLO corrections are in the 
low MET regions, where NLO corrections reach K factors 
of about 1.4 for MET ~ 150 GeV, with a steady decrease 
with increasing MET. We observe similar features also in 
the high-scale benchmark points of 0( 1 ) TeV for my or 
2mx (bottom panels of Fig. [4]), where the largest I\ fac¬ 
tors are about 1.2 for MET ~ 150 GeV. Comparing with 
the FO distributions in Fig. [2j we observe that the parton 
shower does not affect the MET distribution. Note that 
the NLO corrections are different for different MET re¬ 
gions, with the largest NLO corrections occurring in the 
lower MET regions where the rate is the highest. Hence 
the careful estimation of NLO effects is very important for 
accurate LHC studies of DM in each signal region. 

Next, we study the features of jet kinematic distribu¬ 
tions produced in association with DM. Figures !©0 
and [8] show example py and p distributions of the hardest 
and second hardest jets for the four benchmark points as 
in Fig. [4| assuming 

MET >150 GeV. (18) 

Distributions of the hardest jet transverse momentum show 
very interesting features. In Fig. [5] we find that, in all 
benchmark points, the LO distributions match the NLO 
predictions at the peak, i.e. py(j 1 ) ~ 150 GeV, to a very 


good degree. The agreement can be attributed to the im¬ 
posed MET cut in (18), which forces the events into a 
back-to-back configuration of the leading jet and the Y\ 
mediator (on average). We also note that the NLO scale 
uncertainty in the peak region becomes very small com¬ 
pared to the LO estimates. 


The NLO corrections to py(j±) distributions affect not 
only the overall rate, but the shape of the distribution as 
well. In the low -py region, K factors are about 1.2 —1.5. In 
the high py region, we find significant NLO effects again 
for the ( my,mx) = (100,1) and (95, 50) GeV cases (top 
panels), but not for the (my, mx) = (100, 500) and (1000, 
50) GeV cases (bottom panels). We note that the scale un¬ 
certainty does not significantly reduce at NLO in the py 
regions away from the peak, especially for the light medi¬ 
ator and DM case (top panels). Significant differences in 
NLO contributions and theoretical uncertainties in differ¬ 
ent regions of the py(ji) spectrum suggest that the proper 
modelling of the hardest jet differential distributions has 
to go beyond the simple scaling by a constant K factor. 


Apart from the highest py jet which is modelled by the 
hard matrix element, all other jets in the LO simulation 
are generated by the parton shower. By contrast, the NLO 
corrections include real emission diagrams which can con¬ 
tain two hard and well-separated partons in the final state 
as well as virtual corrections to one parton emission. One 
could expect significant differences between LO and NLO 
in the kinematic distributions of the second highest py jet. 
For the (my,mx) = (100,1) and (95,50) GeV cases (top 
panels), we observe giant K factors in the high-p-r tails 
of the distributions. The large difference between LO and 
NLO computations is a consequence of the inadequacy of 
the parton shower to accurately model high-pr emissions. 

In Fig. [6] (bottom panels), on the other hand, we find no 
significant differences between LO and NLO for the over¬ 
all rate and shape of the second jet emission in case of 
very heavy mediators (i.e. my = 1 TeV) or heavy DM 
(i.e. mx = 500 GeV), suggesting that the second hardest 
jet is described very well by the parton shower. This is 
because the scale of the shower is very high and there¬ 
fore extra parton emission from the parton shower can be 
sufficiently hard. 

Features similar to those observed in py(ji^) also oc¬ 
cur in distributions of the hardest/second-hardest jet pseudo¬ 
rapidity (p(ji, 2 )), shown in Figs. [7] and [8] For the light 
mediator/DM (top panels), we observe that the rate at 
which the hardest jet is emitted at NLO in the low ra¬ 
pidity region is enhanced by a factor about 1.5, with the 
corrections falling off with the increase in rapidity. How¬ 
ever, even though the overall rate for the second-hardest 
jet increases by a factor of roughly ~ 1.5 the shape of the 
r/(j 2 ) distribution is affected only mildly. In the case of 
heavy mediator/DM, the hardest jet is emitted at a lower 
rapidity (on average) at a significantly higher rate com¬ 
pared to light mediators as illustrated by the width of the 
p distributions in Fig. [TJ As the hardest jet typically re¬ 
coils against MET, this explains why the MET spectrum 
falls off more quickly for lighter mediators than for the 
heavier ones. 
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Fig. 4. MET distributions for pp — ¥ XX + j at the 13-TeV LHC for four benchmark points specified by (my ,mx), where we 
assume a pure vector mediator and Dirac DM. The middle and bottom panels show the differential scale uncertainties and K 
factors, respectively. 
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Fig. 5. Pt distributions of the hardest jet for pp —> XX + j at the 13-TeV LHC for four benchmark points specified by 
( rriY,mx ), where we assume a pure vector mediator and Dirac DM and the MET > 150 GeV cut is imposed. The middle and 
bottom panels show the differential scale uncertainties and K factors, respectively. 
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Fig. 6. Same as Fig. [5] but for the second hardest jet. 
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Fig. 7. r/ distributions of the hardest jet for pp —> XX +j at the 13-TeV LHC for four benchmark points specified by (my, mx), 
where we assume a pure vector mediator and Dirac DM and the MET > 150 GeV cut is imposed. The middle and bottom 
panels show the differential scale uncertainties and K factors, respectively. 
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Fig. 8. Same as Fig. It| but for the second hardest jet. 
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Fig. 9. Jet multiplicity for pp —> XX+ jets at the 13- 
TeV LHC with FxFx merging in case of my = 1 TeV with 
mx = 50 GeV, where NLO+PS samples are merged up to 
one (red) and two (blue) jets. The 1-jet NLO+PS sample 
without merging (orange) is also shown for comparison. The 
middle panel shows the relative scale uncertainties, while the 
bottom panel presents the ratio of the 2 j merged sample to 
the 1 j merged one and to the non-merged one. We assume 
MET > 100 GeV and jets with pr > 30 GeV and \r/\ < 4.5 for 
the purpose of illustration. 


Fig. 10. Same as Fig. [ 9 ] but for the MET distribution. 

differences are even milder between the samples merged 
to lj and to 2 j. As the 0 j bin is phenomenologically irrel¬ 
evant, we can conclude that the effects of jet merging at 
NLO are within 10%. 

We show effects of jet merging on the MET distri¬ 
bution in Fig. [To] Except in the low MET region, we find 
that the effects of NLO merging are again mild and within 
10 %. 


3.3 Merging samples at NLO accuracy 


In addition to total and differential production cross sec¬ 
tions for the pp -4 XX + j process, we study NLO ef¬ 
fects for different jet multiplicities in the final state. For 
this purpose we utilise the FxFx merging procedure 
within the framework of MG5 aMC, and consider 


20 


pp -+ XX + 0,1,2 jets. 


(19) 


We take the merging scales at 45 and 30 GeV for 2- and 
1-jet merged samples, respectively. 

Figure [9] shows the number of jets in the final state 
for NLO merged samples in case of my = 1 TeV and 
m.\ = 50 GeV. The red and blue curves show the re¬ 
sults of merging up to 1 and 2 jets, respectively, while 
the orange curve shows the pp -+ XX + lj process at 
NLO+PS without merging for comparison. An inspection 
of the three samples in the lowest panel of Fig. [9] shows 
that the effects of NLO merging are mild. The non-merged 
NLO sample over-estimates the production rate in the 0 j 
and > 2 j bins by 20 and 10% respectively, and under¬ 
estimates the rate in the 1 j — 2 j bins by < 10 %. The 


3.4 Comparison of signal distributions to the Standard 
Model 

In discussions of NLO corrections to DM production, it 
is important to consider how the possible signal events 
at NLO look in the midst of large SM backgrounds. For 
the purpose of illustration, we consider only the largest 
background in mono-jet searches for DM, i.e. Z +jets and 
simulate it to NLO merged up to 2 extra jets via the FxFx 
method. 

Figure [Tl] shows an example comparison of the Z+jets 
channel to several benchmark points of the simplified model 
discussed in previous sections. The shape of signal jet 
multiplicity distributions (upper left panel) resembles the 
Z +jets distribution to a good degree, while the overall 
rate varies wildly depending on the model point. Note 
that events containing one jet are produced at almost an 
identical rate to 2 j events and comparable to 3 j events, 
both in Z +jets and all of the benchmark model points 
we considered. The production rate for different jet multi¬ 
plicities implies that it could be beneficial to consider DM 
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Fig. 11. Distributions for various signal benchmark points and Z(uu)+ jets, where MET > 150 GeV and jets with pr > 30 GeV 
and I 77 I < 4.5 are considered. We assume a pure vector mediator and Dirac DM with the coupling parameters gx = 1 and 
gSM = 0.25. 


searches beyond mono-jet, either inclusive, or at fixed jet 
multiplicities. 


4 Dark matter production with a top-quark 
pair 


Next, the MET (as well as the hardest jet pp) distri¬ 
butions for mediators of mass ~ 100 GeV with light DM 
naturally resemble the Z +jets distribution in shape, as the 
kinematics of the MET and hardest jet are determined by 
the mass scale of the heavy object ( i.e. the Z boson or the 
mediator). For heavy mediators, the MET andp-rC/i) dis¬ 
tributions fall off with a milder slope, suggesting that the 
signal to background ratio ( S/B ) in searches for DM could 
be improved by requiring a higher for the hardest jet. 

The pt distribution of the second hardest jet, on the 
other hand, seems to display a similar shape in all model 
points, as well as the Z +jets background channel. Requir¬ 
ing a high pt on the second jet would hence not improve 
neither S/B nor the signal significance, suggesting that in 
searches which exploit the presence of a second jet, only 
minimum cuts on pr(j 2 ) should be applied. Note that the 
modelling of the second jet from LO calculations grossly 
underestimates the overall rate in case of lighter mediators 
f° r Pt > 100 GeV, as illustrated in Fig. [bj 


In the spin-0 mediator model, due to the normalisation 
of the Yukawa couplings in the Lagrangian (111, the most 
relevant tree-level process at the LHC is 


pp —> XX + tt. 


( 20 ) 


Such models have in the past been studied in the con¬ 
text of EFT 157-59 and simplified models 146 48 , and 


searched for at the LHC Run I 60,61 . Past work on DM 


interactions with the top quarks has mainly focused on LO 
estimates, with only a few analyses including NLO correc¬ 
tions. Here we present a comprehensive study of NLO ef¬ 
fects of DM interaction with top quarks in the framework 
of the simplified model. We note that a wide class of so 
called “top-philic DM” models exist where the LO pro¬ 
duction is via top loops. Reference 62 studied such a sce¬ 


nario in a minimally model-dependent framework, while 
more recently, Ref. [ 44] presented concrete predictions for 
loop-induced DM production for the current LHC13 run, 
using the same simulation framework as in this work. 
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(■ m Y ,m x ) [GeV] 



scalar 

pseudo-scalar 

10 

undecayed 

0 lo [pb] 
onlo [pb] 
K factor 

2.278 x 10 1 ^20.4±4-2% 
2.435 x lO^g.tii.sro 

1.07 

5.202 x 10- 1 lf 2 :o±6.o% 
5.431 x 10~ 1 tlo 4 2±2.6% 
1.04 

(10, 1) 

my > 2 m x 

OLO [pb] 
onlo [pb] 
K factor 

2.294 x 10 1 ^2o.5±4-2% 
2.460 x 10 1 ±1.8% 

1.07 

5.500 x HP 1 ± 6 . 0 % 

5.739 x 10' 1 llo 4 2±2.6% 
1.04 

(10, 50) 

my < 2 m x 

olo [pb] 
ONLO [pb] 
K factor 

2.415 x 10 —3 — 2 i".s ± 5.8 % 
2.340 x 10 " 3 ^ 9 .?± 2 . 8 % 
0.97 

3.329 x 10 " 3 tlH ± 8 . 7 % 
3.133 x 10 -3 1 ilo ±3.9 % 
0.94 

100 

undecayed 

OLO [pb] 

onlo [pb] 
K factor 

8.226 x 10 _1 — 20.9 ± 4 . 4 % 
8.391 X 10 — 1 1^8.6 ±2.1 % 
1.02 

2.442 x 10 ” 1 ± 22(9 ±7.2% 
2.431 x 10 — 1 —icl7 ±3.2 % 
1.00 

( 100 , 1 ) 

my> 2 m x 

olo [pb] 

ONLO [pb] 
K factor 

8.135 x 10 _1 — 20.9 ±4.4% 
8.207 x 10 _ 1 — 8.3 ± 2.1 % 
1.01 

2.464 x 10” 1 ±7.2% 

2.427 x 10~ 1 llo°4±3.2% 
0.98 

(95, 50) 

my < 2 m x 

O' LO [pb] 
ONLO [pb] 
K factor 

7.986 x 10 " 3 ^ 2 i'. 3 ± 5 .o% 
7.897 x 10 —3 1^8.8 ± 2.4 % 
0.99 

1.404 x 10 " 2 ± 23(3 ± 7 . 8 % 
1.362 x 10 -2 1I 0 4 8 ± 3 . 5 % 
0.97 

1000 

undecayed 

OLO [pb] 

onlo [pb] 
K factor 

1.571 x 10 —3 — 2 t". o ±17.1 % 
1.127 X 10 " 3 ^3.6 ±8.7% 
0.72 

1.827 x 10 " 3 ± 2 ?;t±i 7 . 4 % 
1.297 x 10 " 3 ll 3 ;?± 8 . 8 % 
0.71 

( 1000 , 1 ) 

my > 2 m x 

olo [pb] 

ONLO [pb] 
K factor 

7.499 x 10~ 4 ^27.2±i6.8% 
5.201 x 10 " 4 ^i 2 4 7 ± 8 . 4 % 
0.69 

8.174 x 10- 4 ^±17.0% 
5.675 x 10" 4 li2 6 9 ±8.5% 
0.69 

(1000, 50) 

my > 2 m x 

O' LO [pb] 

onlo [pb] 
K factor 

7.354 x 10 _ 4 ^ 27.2 ± 16.8 % 
5.125 x 10" 4 ^2 6 s±8.5% 
0.69 

8.137 x 10 -4 ^ 27(3 ±17.0% 
5.595 x 10 - 4 ti 2 3 7 ± 8 . 5 % 
0.69 


Table 3. LO and NLO cross sections and corresponding K factors for DM pair production in association with a top-quark 
pair for the scalar and pseudo-scalar mediator scenario at the 13-TeV LHC. The uncertainties represent the scale and PDF 
uncertainties in per cent, respectively. We show several benchmark model points for the mediator and DM masses with the 
coupling parameters g x = 1 and gsM = 1. 


The code and events for the above process can be au¬ 
tomatically generated by issuing the following commands 
in MG5AMC: 


> import model DMsimp_s_spinO 

> generate p p > xd xd~ t t~ [QCD] 

> output 

> launch 


We have checked that our model can reproduce the SM 
predictions for pp -t htt —» T + r~tt by adjusting the cou¬ 
pling and mass parameters. Note that we use the on-shell 
renormalisation for the NLO model construction. The top- 
quark decays can be subsequently performed by Mad- 
Spin 63 , which keeps production and decay spin correla¬ 
tions. 

To illustrate the NLO effects, we consider pure scalar, 
Eqs. (12) and (13), or pure pseudo-scalar, Eqs. (141 and 
(15), couplings, and take 


( 21 ) 


as the default couplings. With these values, the scalar and 
pseudo-scalar mediator width is ry/my ~ 0.06 — 0.1 for 
my > 2mx,2mt, while fy/my ~ 0.04 for my > 2m x 
and my- < 2 m t - The width for scalar is slightly smaller due 
to the additional / 3 2 factor where /3 = (1 —Am\ t /my) 1 ^ 2 . 


4.1 Total cross sections 

We start by showing total production rates of pp —> XX + 
ti, where the top quark is considered stable. Table[3]shows 
the LO and NLO cross sections (in pb) for the scalar 
and pseudo-scalar mediator scenarios, where we use mt = 
172 GeV. The central renormalisation and factorisation 
scales are set to half the sum of the transverse mass of 
the top quarks and the missing transverse energy. We also 
present scale and PDF uncertainties in % as well as K 
factors. 

For the total production rates, the NLO effects are very 
mild for the light mediator case, while they are significant 


gx = 1 and g SM = 1 
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Fig. 12. Distribution of the invariant mass of the top-quark pair for pp —> XX + tt at the 13-TeV LHC for different mediator 
masses with the DM mass fixed at 50 GeV, where we assume a pure scalar mediator and Dirac DM. The middle and bottom 
panels show the differential scale uncertainties and K factors, respectively. 





300 400 500 600 700 800 


300 400 500 600 700 800 


300 400 500 600 700 800 


m(tt) [GeV] 


m(tt) [GeV] 


m(tt) [GeV] 


Fig. 13. Same as Fig. 


12 


but for the pseudo-scalar mediator scenario. 


for the heavy case. The inclusion of NLO corrections re¬ 
sults in a drastic improvement of the scale uncertainties, 
from up to 40 % at LO to typically only about 10 % at 
NLO. Also, the PDF uncertainties are reduced by approx¬ 
imately a factor of two when going from LO to NLO. 

Table! also shows clear differences between the over¬ 
all production rates in the cases of scalar and pseudo¬ 
scalar mediators. For mediator mass of 0(10) GeV we 
find that DM production cross section via scalar medi¬ 
ators is an order of magnitude larger compared to the 
production rate via the pseudo-scalar mediator with the 
same mass. The large difference occurs due to the fact 
that in case of rrix < wy <C mt, the production cross 
section is dominated by the t —> tYo fragmentation. In 
case of the scalar mediator, the t —► tYo fragmentation 
function contains terms with soft singularities of the form 
(1 — x)/x - where x is the momentum fraction carried by 
the mediator - causing enhancements in the production 


rate 64 . The soft-enhanced term is absent in the case 


of a pseudo-scalar mediator 65 , explaining the order of 
magnitude difference between the total rates of the scalar 
and pseudo-scalar mediators. 


In cases where either DM or the mediator is produced 
close to threshold, we observe that the production cross 
section in the pseudo-scalar mediator case is larger. The 
effect can be attributed to the production rate originating 
mainly from top fusion diagrams. The production of a 
DM (Dirac) pair via scalar mediators tt Yq XX 
at threshold can proceed only via a P -wave ( 3 P 0 ) and is 
hence suppressed by extra two powers of /I = \J 1 — 4m^ /s 
661. Conversely, production of DM pair via pseudo-scalar 
mediators can proceed via an S -wave ( 1 S'o) and hence does 
not suffer any kinematic suppression. 
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Fig. 14. Comparison of tt invariant mass distributions between 
the scalar (blue) and pseudo-scalar (red) mediator models for 
different mediator masses. 


4.2 Differential distributions 


For the study of differential distributions, we consider the 
invariant mass of the top-quark pair ( m(ti )), without in¬ 
clusion of a parton shower. Figures p~2| and p~3| illustrate the 
scalar and pseudo-scalar results, respectively, for different 
mediator masses (off-shell, threshold, and on-shell) with 
the DM mass fixed at 50 GeV. In all cases the shape of 
m(tt) is well modelled by the LO calculation, and includ¬ 
ing a constant K factor reproduces the NLO results to an 
excellent degree, except in the threshold region. However, 
the LO calculation suffers from significant scale uncertain¬ 
ties which tend to increase with m(tt), whereas the scale 
uncertainties are under much better control at NLO. 

Whether DM is produced via scalar or pseudo-scalar 
mediators can have a dram atic effect on the shape of the 
m(tt) distribution. In Fig. 14 we compare the NLO dis¬ 


tributions in Figs. |12| and [l3[ where we normalise the his¬ 
tograms to unit area to point out the shape differences. 
We observe that the shape of the distribution is partic¬ 
ularly enhanced for rn(ti) > 500 GeV in the case of the 
pseudo-scalar mediator, while the scalar mediator distri¬ 
bution displays a much more prominent peak at lower 
m(tt). The effect is severely damped in case of heavy me¬ 
diators, where we find no clear differences between the 
shapes of the scalar and pseudo-scalar mediator distribu¬ 


tions. Figure [L4] suggests that scalar and pseudo-scalar 
mediators could be distinguished based on the shape of 
the m(tt) distribution, as long as the mediator is suffi¬ 
ciently light and/or does not decay highly off-shell. An 
analogous observation has been made already in the case 
of the study of the CP properties of the Higgs boson 67 . 


5 Summary 

Searches for DM are one of the main endeavours at the 
LHC Run II. Accurate and precise predictions for produc¬ 
tion rates and distributions are necessary to obtain robust 


constraints on DM models and characterise possible DM 
signals. In this article we have provided a general imple¬ 
mentation of the simplified DM model approach into a 
calculation/simulation framework that allows to system¬ 
atically evaluate and include NLO QCD corrections to 
the production of DM at the LHC. We have considered 
a class of simplified models where DM is a Dirac fermion 
and couples to the SM via either spin-1 or spin-0 s-channel 
mediators, making no restrictions on chiral couplings. For 
the purpose of illustration, we analysed the NLO effects on 
the DM production via vector and axial-vector mediators 
in the context of mono-jet signals. In addition, we have 
presented detailed predictions of DM production in asso¬ 
ciation with a top-quark pair via scalar and pseudo-scalar 
mediators. We presented our results for various DM and 
mediator masses to cover benchmark points suggested by 
the ATLAS/CMS DM forum [5]. 

For MET+jets in the spin-1 mediator model, our re¬ 
sults show that higher-order corrections have a significant 
effect both on the overall production rate as well as on the 
shape of relevant differential distributions, with a sizeable 
reduction of the scale and PDF uncertainties. The NLO 
corrections to the LO production rates can be large, with 
K factors of up to K < 2, and typically occur in parts of 
the model parameter space where the mass scale of DM 
and mediator is 0(10 — 100) GeV. For such scenarios, 
we also find large NLO effects on the shape of differen¬ 
tial distributions in MET and the transverse momentum 
of the associated jets. Simplified models with heavy ( e.g. 
0(1) TeV) mediators/DM do not receive large NLO cor¬ 
rections, and we find that LO predictions describe both 
total production rates and shapes of differential distribu¬ 
tions quite accurately. Distributions of the second hardest 
jet in the event are well modelled by the parton shower 
for heavy mediator/DM cases. On the other hand, for me¬ 
diators/DM with masses of 0(100) GeV, the inclusion of 
NLO effects is essential for a proper description of pr(j 2 ) 
and i)(j 2 ) distributions, especially in the high-py tails, 
where the NLO effects can be an order of magnitude. 

So-called “giant K factors” can occur in NLO com¬ 
putations of DM production rates in the regions where 
p J T my. Such effects can be extremely large when con¬ 
sidering mono-jet production rates, especially in phase- 
space regions with low MET. Imposing a sufficiently large 
MET cut and hence avoiding the soft/collinear singulari¬ 
ties associated with the mediator emissions from high -px 
jets efficiently mitigates the effect of giant K factors. 

In our analysis we have gone beyond FO in pertur¬ 
bation theory and studied the effects of jet multiplicity 
merging at NLO accuracy. We found that FO calculations 
model the jet multiplicity and other differential distribu¬ 
tions adequately well, with no significant effects on the 
shapes or overall rates coming from jet sample merging. 

Comparisons with the NLO predictions for the leading 
SM background channel reveal that considerations of ei¬ 
ther inclusive or exclusive jet samples beyond one jet could 
be beneficial for increasing the prospects for DM detec¬ 
tion. The leading jet pt distributions in case of heavy me¬ 
diators display a milder decrease with the increase in pr, 
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suggesting that a significant improvement in S/B could 
be obtained by focusing on high-p^ regions. The second 
hardest jet pt distribution, on the other hand, was char¬ 
acterised by a slope similar to the leading backgrounds, 
implying that more inclusive cuts on the second jet should 
be used. 

For MET+ft in the spin-0 mediator model, our re¬ 
sults show that the NLO corrections are very mild for 
the light mediator case, while they are significant for the 
heavy case. We observed a drastic improvement of the the¬ 
ory uncertainties when going from LO to NLO. We have 
noted that the shape of the m(tt) distribution can reveal 
the chiral structure of the DM-SM interactions, as long 
as the mediator is relatively light ( i.e. < (9(100) GeV). 

The DM model we studied in this paper is publicly 
available in the FeynRules repository [28], We emphasise 
that all results presented here have been obtained in the 
FeynRules/MG5aMC framework, and thus they can be 
easily reproduced and used in DM searches at the LHC 
Run II. 
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